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ABSTRACT 

Surface plasmon resonance was used to investigate two previously described interactions analyzed by reverse genetics and 
complementation mutation experiments, involving 5BSL3.2, a stem-loop located in the NS5B coding region of HCV. 5BSL3.2 
was immobilized on a sensor chip by streptavidin-biotin coupling, and its interaction either with the SL2 stem-loop of the 3' 
end or with an upstream sequence centered on nucleotide 9110 (referred to as Seq9110) was monitored in real-time. In 
contrast with previous results obtained by NMR assays with the same short RNA sequences that we used or SHAPE analysis 
with longer RNAs, we demonstrate that recognition between 5BSL3.2 and SL2 can occur in solution through a kissing-loop 
interaction. We show that recognition between Seq9110 and the internal loop of 5BSL3.2 does not prevent binding of SL2 on 
the apical loop of 5BSL3.2 and does not influence the rate constants of the SL2-5BSL3.2 complex. Therefore, the two binding 
sites of 5BSL3.2, the apical and internal loops, are structurally independent and both interactions can coexist. We finally show 
that the stem-loop SL2 is a highly dynamic RNA motif that fluctuates between at least two conformations: One is able to 
hybridize with 5BSL3.2 through loop-loop interaction, and the other one is capable of self-associating in the absence of 
protein, reinforcing the hypothesis of SL2 being a dimerization sequence. This result suggests also that the conformational 
dynamics of SL2 could play a crucial role for controlling the destiny of the genomic RNA. 
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INTRODUCTION 

The Hepatitis C virus (HCV) is a human viral pathogen of the 
family of the Flaviridae. The HCV genome is a positive-sense 
single-stranded RNA that encodes a long open reading frame 
(ORF). This ORF is translated as a polyprotein of about 3010 
amino acids (Clarke 1997). The coding sequence is flanked 
by two untranslated regions (UTR). The 341-nt 5' UTR and 
the first nucleotides encoding the capsid protein form a highly 
structured domain that serves as an internal ribosome entry 
site (IRES) initiating the ORF translation in a cap-indepen- 
dent manner. The 3' UTR is composed of a short variable re- 
gion, a polyU/C tract of variable length and a 98-nt sequence 
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(3'X), highly conserved among various isolates (Tanaka et al. 
1995), which folds as three contiguous stem-loops, SLl, SL2, 
and SL3 (Fig. lA; Friebe and Bartenschlager 2002; Yi and 
Lemon 2003). 

Evidence has now accumulated suggesting that RNA-RNA 
interactions between structured regions of the HCV genome 
are crucial to regulate the life cycle of the virus (You et al. 
2004; Friebe et al. 2005; Diviney et al. 2008; You and Rice 
2008; Tuplin et al. 2012). 

Several interactions in the HCV genome involve a struc- 
tured RNA motif, known as 5BSL3.2, located within the 
NS5B coding region (Fig. lA). This imperfect hairpin of 48 
nt has been shown to play a role in the life cycle of HCV and 
would act as a molecular hub capable of interacting with other 
RNA regions within the viral RNA. 5BSL3.2 displays a 7-nt 
sequence in an apical loop (Fig. IB) that would recognize 
stem-loop SL2, through kissing-loop interactions (You et al. 
2004; Friebe et al. 2005; You and Rice 2008). Molecular re- 
cognition between these two RNA modules has been shown 
to be critical for replication of the replicon (You et al. 2004; 
Friebe et al. 2005). The internal loop of 5BSL3.2 contains 
6 nt which would hybridize with a sequence centered on 
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FIGURE 1. (A) Scheme of the HCV genome. The internal ribosome 
entry site (IRES) is on the 5' side. The SL2 stem-loop is located in the 
3' UTR region. Seq9110 and 5BSL3.2 are located in the region coding 
for the NS5B protein. The region coding for the core protein (Core) 
is also shown. (B) Sequences and secondary structures of the RNAs. 
5BSL3.2 was modified with a biotin tag at the 3' end, allowing immobi- 
lization by streptavidin-biotin coupling on the sensor chips. The se- 
quence in bold black indicates bases that could form Watson-Crick 
base pairs with SL2 (also in bold black). The sequence in bold gray in- 
dicates bases that could interact with Seq9110 or ASgno (also in bold 
gray). Mini5BSL3.2 is a truncated version (double arrow) of 5BSL3.2 
corresponding to the perfect stem-loop above the internal loop. 
AS9110 is an antisense oligonucleotide derivatized from Seq9110. For 
SL2, the sequences in italic gray indicate the palindromic sequence of 
SL2. Underlined bases for SL2 and Seq9110 indicate nucleotides that 
were replaced to generate hairpins of decreased affinity for 5BSL3.2: 
UGU was replaced by ACA for SL2, and GGG by AAA for Seq9110. 
SL2wTi SL2wT2 are two predicted secondary structures of SL2 ob- 
tained from the Mfold web server (http://mfold.rna.albany.edu/? 
q=mfold/RNA-Folding-Form). Bases in bold gray for SL2]jjs(en, (A and 
C) indicate modifications that were made to convert SL2wt2 in a perfect 
stem-loop structure. SL2gc was designed to generate a hairpin of in- 
creased stability compared with SL2 without altering the palindromic 
sequence. SL2aau is a mutated version of SL2gc in which CUA was mu- 
tated in AAU (in bold gray). 



nucleotide 9110 (referred to as Seq9110) (Fig. IB; Diviney 
et al. 2008). This internal loop has also been implicated in a 
long-range RNA-RNA interaction with stem-loop Illd of 
the IRES (Romero-Lopez and Berzal-Herranz 2009), but 
there are conflicting results about the influence of this interac- 
tion on viral translation activity. It is hypothesized that 
5BSL3.2 could act as a molecular switch depending on its in- 



teraction with SL2 or Seq91 10, therefore leading the viral ge- 
nome toward translation, replication, or encapsidation 
(Tuplin et al. 2012). 

Despite the role that these interactions are expected to play 
in the life cycle of the virus, they have not as yet been charac- 
terized by direct kinetic methods. Some of them could not 
even be observed when they were analyzed by NMR spectro- 
scopy. We, therefore, wanted to address whether 5BSL3.2 and 
SL2 could, indeed, form a loop-loop complex, whether 
Seq9110 could interact with the internal loop of 5BSL3.2, 
and whether these interactions were mutually exclusive or 
not. For this purpose, we used surface plasmon resonance 
(SPR). This label-free technique can be used to assess the af- 
finity and the kinetics of bimolecular (Morton and Myszka 
1998) or even trimolecular complexes (Bernat et al. 2003; 
Walsh et al. 2003; Jomain et al. 2007; Di Primo 2008; Palau 
and Di Primo 2012). In some cases, stoichiometry can also 
be determined (Teulade-Fichou et al. 2003; Mistrik et al. 
2004; Di Primo and Lebars 2007; Di Primo 2008). SPR is 
widely used to investigate protein-protein, protein-small mol- 
ecules, or protein-DNA interactions (for review, see Rich and 
Myszka 2011) but rarely for investigating RNA-RNA com- 
plexes, likely due to their susceptibility to hydrolysis. In pre- 
vious studies, we have successfully investigated interactions 
between RNA oligonucleotides with structured RNA of hu- 
man viruses or cells (Duconge et al. 2000; Aldaz-CarroU et al. 
2002; Lebars et al. 2008; Dausse et al. 2011). 

In this work, the imperfect stem-loop structure 5BSL3.2 
was immobilized on a streptavidin-coated sensor chip. SL2 
or Seq9110 RNA oligonucleotides were injected in a contin- 
uous flow of buffer. These interactions were characterized 
either in a bimolecular context with the 5BSL3.2-SL2 and 
5BSL3.2-Seq9110 complexes or in a trimolecular context 
with the 5BSL3.2-SL2-Seq9110 ternary complex. In contrast 
with previous results obtained by NMR spectroscopy with 
the same short RNA sequences that we used in this work 
(Friebe et al. 2005) or SHAPE mapping with longer RNAs 
(Tuplin et al. 2012), we show that 5BSL3.2 can form a com- 
plex with SL2 through kissing- loop interactions. The internal 
loop of 5BSL3.2 is also able to hybridize with Seq9110. The 
kinetic analysis of the 5BSL3.2-SL2-Seq91 10 ternary complex 
shows that the two bimolecular complexes 5BSL3.2-SL2 and 
5BSL3.2-Seq9110 are not mutually exclusive. The two inter- 
actions can coexist, strongly suggesting that the two binding 
sites of 5BSL3.2, the apical and the internal loop, are structur- 
ally independent. We also found that SPR signals that could 
not be fitted with a 1:1 model of interaction when SL2 was 
injected over 5BSL3.2 could only be explained by the struc- 
tural heterogeneity in a solution of SL2. We demonstrate 
that this stem-loop is, indeed, a highly dynamic structure 
that alternates between at least two conformations: one is 
able to recognize 5BSL3.2 through kissing interactions, and 
the other one is able to self-associate. Both of these interac- 
tions occur in the absence of any added protein. This result 
suggests that the conformational dynamics of SL2 could 
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play a crucial role for controlling the destiny of the genomic 
RNA. 



RESULTS 

Design of the RNA molecules 

The RNA molecules used for the SPR experiments were de- 
signed according to previously published work on HCV ge- 
notype lb (isolate Conl) (You et al. 2004; Friebe et al. 
2005; You and Rice 2008). Located within the NS5B coding 
region, 5BSL3.2 folds as an imperfect hairpin displaying an 
apical and an internal loop (Fig. IB) and has been shown 
to be crucial for replication. The predicted structure by the 
mfold server (Zuker 2003) was confirmed by NMR spectro- 
scopy (Friebe et al. 2005). 5BSL3.2 was modified with abiotin 
tag at its 3' end to allow immobilization on a streptavidin- 
coated sensor chip. SL2 was also designed according to the 
known secondary structure determined within the 3'X of 
HCV (Friebe and Bartenschlager 2002; Yi and Lemon 
2003). This stem-loop was referred to as SL2 wild type 
(SL2wt)- The mfold server predicted three secondary struc- 
tures for SL2. Only those that displayed the sequence for rec- 
ognizing 5BSL3.2, partially (SL2wti> AGsy-c = — 8.20 kcal/ 
mol) or entirely (SL2wt2> AGsy-c = -7.60 kcal/mol) in an api- 
cal loop are reported in Figure IB. SL2 has been previously 
characterized by NMR spectroscopy (Friebe et al. 2005). In 
this work, the authors used 1.4 mM of RNA prepared either 
in 20 mM NaCl-5 mM acetic buffer, pH 5.5, or in 20 mM 
NaCl-20 mM phosphate buffer, pH 7. In these experimental 
conditions, only conformation SL2wt2 was observed in solu- 



tion. Seq9110 (Fig. IB) was designed as a 2 1-nt-long RNA dis- 
playing the sequence that potentially recognizes the internal 
loop of 5BSL3.2 at a central position. The mfold server pre- 
dicted one secondary structure (Fig. IB) with an initial 
AG37-C close to zero (—0.9 kcal/mol). The region encompass- 
ing nucleotide 9 1 1 0 is supposedly unstructured (Diviney et al. 
2008). To check whether or not a structured context was cru- 
cial for recognition, an antisense oligonucleotide was also de- 
signed, ASgiio (Fig. IB). This oligonucleotide contained the 
6-nt recognition sequence and was flanked by CA nucleotides 
that are not expected to contribute to binding, in order to in- 
crease the SPR signal by increasing its molecular weight. As 
expected, no structure could be predicted by mfold. 

Binding of SL2 and Seq91 1 0 to 5BSL3.2 

The results obtained when SL2wt was injected over the 
5BSL3.2 functionalized surface are reported in Figure 2A. 
These experiments were performed using the SCK method, 
as described in Materials and Methods, which is faster than 
the classical method and avoids repeated steps of injection/ 
regeneration that could alter the surface and the immobilized 
target. Binding is observed when SL2wt is injected at increas- 
ing concentrations. Flat sensorgrams were obtained when a 
mutated stem-loop in which UGU (underlined in Fig. IB) 
was replaced by ACA was injected {SLl^^J (Fig. 2A). This 
strongly confirms that the interaction between 5BSL3.2 and 
can occur through kissing interactions between the se- 
quences that have been previously proposed (Lee et al. 2004; 
You et al. 2004; Friebe et al. 2005). Interestingly, the sensor- 
grams (Fig. 2 A) are not monophasic. When the results are 
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FIGURE 2. Kinetic analysis by the SCK method of SL2 hairpins binding to 5BSL3.2. The experiments were performed at 10°C on a SAHC200m 
sensor chip, at a flow rate of 50 j^L/min. Thirty to 70 RU of biotinylated 5BSL3.2 were immobilized on the sensor chip. The samples were injected 
sequentially in the order of increasing concentrations, 62.5 nM (first arrow from the left), 250 nM (second arrow), and 1000 nM (third arrow). 
The regeneration was achieved with a 2-min pulse of a mixture of 40% formamide, 3.6 M urea, and 30 mM EDTA prepared in milli-Q water. 
The gray curves represent the experimental data (two overlaid sensorgrams). The black line represents the fit of one sensorgram to a Langmuir 
1:1 model according to Equations (1) and (2) reported in the Supplemental Material. (A) SL2 wild type (SL2wt) and the mutated SL2 stem-loop 
(SL2aca) in which the central UGU sequence of the loop was replaced by ACA. (B) SL2]5is,en,. 
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fitted with a simple Langmuir 1:1 model, there is a clear dis- 
crepanq^ between the data and the model used to fit them. 
Instead of trying to analyze the data with a more complicated 
model assuming, for instance, a conformational change after 
binding, a popular model often used to fit complex signals, 
we hypothesized that the nonmonophasic shape of the sen- 
sorgrams was simply due to a structural heterogeneity of 
SL2\YT- 

We then analyzed SL2wt RNA using native gel electropho- 
resis in buffer conditions similar to those used for the SPR ex- 
periments, in particular, in the presence of 3 mM magnesium 
cations, in order to determine whether several structures do 
exist in solution (alternative monomeric structures, dimers) 
(Fig. 3). Strikingly, in contrast with previous results obtained 
by NMR spectroscopy suggesting that SL2 adopts only one 
conformation in solution, several bands and a smear are ob- 
served for SL2wT) in the presence of Mg^^ and at much lower 
concentration of RNA (45 |-iM for the loaded sample instead 
of >1 mM for the NMR assays). This indicates that SL2wt 
may rapidly alternate between several structures, monomers 
(lower bands), and dimers (upper bands). The hypothesis 
that the unusual sensorgrams could be due to the hetero- 
geneity of the 5BSL3.2 oligonucleotide was discarded be- 
cause a single migrating band was visualized after loading 
of 5BSL3.2 on a native acrylamide gel (Supplemental Fig. SI). 

We then hypothesized that mutations of SL2wt that would 
generate a perfect hairpin displaying in its apical loop the 
sequence for recognizing 5BSL3.2 by loop-loop interactions 
should result in one single band on the native gel, and the 
sensorgrams might be fitted with a 1 : 1 model. For this pur- 
pose, two modifications were made on the SL2wt2 confor- 
mation (Fig. IB). First, an adenine was added to generate 
an AU Watson- Crick pair with the unpaired U in the wild- 
type conformation (SL2wt2)) and the UG pair in the stem 
next to the loop was changed to CG. The mfold web server 
predicted one unique stem-loop structure (AGsy-c = —14.80 
kcal/mol), referred to as SL2kistem in Figure IB. One single mi- 

-V*- ^v* 'V^ 'V'^ 

cs' or 
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FIGURE 3. Native acrylamide gels of SL2 hairpins. Two of SL2wt, 
SL2ids,eni, SL2gc! 3nd SL2aau were prepared and were loaded on 15% 
(w/v) 75:1 acrylamide/bis(acrylamide) native gels as described in 
Materials and Methods. The gels were stained by "Stains-all." 



grating band was observed when SL2kis,j.m was loaded on a 
native gel (Fig. 3), and the sensorgrams were perfectly fitted 
with a Langmuir 1:1 model (Fig. 2B). The rate and the dis- 
sociation equilibrium constants for the 5BSL3.2-SL2ystem 
kissing complex are reported in Supplemental Table SI. 
In this work, only SL2i;is,j.rn will be used for further analyses 
with 5BSL3.2. 

We then investigated the interaction between 5BSL3.2 and 
Seq9110. In the preliminary experiments, RNA-RNA bind- 
ing was characterized using the SCK method (Supplemental 
Fig. S2). Clear interaction was observed with slow association 
and dissociation phases. However, attempts to fit the sensor- 
grams with a Langmuir 1 : 1 model were not conclusive. The 
curve fittings clearly deviated from the sensorgrams, in par- 
ticular, for the association phases. The experiments were re- 
peated using the classical SPR method to analyze carefully the 
association phases on longer time scales (Fig. 4A). In these 
experiments, three concentrations of Seq9110 were injected 
in duplicate. The regeneration step was performed at least 
800 sec after the beginning of the dissociation phase. Either 
at low, 250 nM, or high concentrations, 5 [iM, duplicate over- 
laid sensorgrams demonstrated that regeneration was effi- 
cient. Again, the data still did not fit to a 1:1 model. Other 
attempts consisting of lowering the level of immobilization 
of biotinylated 5BSL3.2 to exclude problems related to a 
too high density of the immobilized target did not improve 
the curve fitting (data not shown). Seq9110 was loaded on 
a native gel to test a potential structural heterogeneity respon- 
sible for these findings. Seq9110 is slightly heterogeneous 
compared to the immobilized target 5BSL3.2 (Supplemental 
Fig. SI) but to a lesser extent than SL2wt (Fig- 3). We actually 
cannot conclude if this is a reason good enough to explain 
why these data could not be fitted with a 1 : 1 model. Despite 
the problems in determining the corresponding rate con- 
stants, the observed kinetics are clearly slower than those 
obtained with SL2ijstem- In addition, the results show that 
the interaction between the internal loop of 5BSL3.2 and 
Seq91 10 is specific. Indeed, substitution of the three guanines 
by three adenines (Seq9110aaa) in the recognition sequence 
of Seq9110 (Fig. IB) drastically decreases the binding (Sup- 
plemental Fig. S3). 

An antisense oligonucleotide, ASgno (Fig. IB), was then 
injected over the immobilized 5BSL3.2 stem-loop (Fig. 4B) 
in order to investigate whether the 6-nt recognition sequence 
of Seq91 10 required a secondary structure context to hybrid- 
ize with the internal loop of 5BSL3.2. Clear binding was 
observed, but again, the data could not be fitted with a 1:1 
model of interaction. 

Since 5BSL3.2 has two binding sites for recognizing other 
RNA sequences of the HCV genome, the next step was to 
investigate the effect of the 5BSL3.2-SL2 interaction on 
5BSL3.2-Seq9110 binding and whether these interactions 
were simultaneous or mutually exclusive. The experiments 
were performed by the single cycle kinetics on a decaying 
surface (SCKODS) method that was developed to analyze 
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FIGURE 4. Kinetic analysis of Seq9110 and ASgno binding to 5BSL3.2. The sensor chip was prepared as described in Figure 2 and in the Materials and 
Methods. RNA samples were injected over the functionalized surface in duplicate. Two overlaid sensorgrams (gray curves) are shown for each injected 
concentration. The regeneration of the surface was achieved with a 2-min pulse of a mixture of 40% formamide, 3.6 M urea, and 30 mM EDTA pre- 
pared in milli-Q water. The black lines represent the fit of one set of sensorgrams by global fitting analysis to a Langmuir 1:1 model. (A) Injection of 
Seq9110. (B) Injection of AS9110. 



ternary complexes (Fig. 5; Palau and Di Primo 2012). 
Seq9110 was first injected over the biotinylated 5BSL3.2 tar- 
get, immobilized on the streptavidin sensor chip. Then, while 
Seq9110 slowly dissociated from 5BSL3.2, the perfect stem- 
loop SL2kistem was injected sequentially in the order of in- 
creasing concentrations (Fig. 5). The sensorgrams are similar 
to those obtained when SL2]jistem was injected on 5BSL3.2 
in a bimolecular context, i.e., in the absence of Seq9110 
(Fig. 2B). The rate and dissociation equilibrium constants de- 
termined by direct fitting of the sensorgrams are similar to 
those obtained when SL2]jis,eni was injected in the absence 
of Seq9110 (Supplemental Table SI). No binding was ob- 
served when SL2aca was injected (Supplemental Fig. S4), in 
agreement with previous results (Fig. 2A). The results indi- 
cate that the two binding sites of 5BSL3.2 (apical and the in- 
ternal loops) are structurally independent of each other. The 
reverse experiment confirms this finding. Saturation of the 
apical loop of 5BSL3.2 by SL2 did not prevent recognition be- 
tween Seq91 10 and the internal loop of the immobilized tar- 
get (Supplemental Fig. S5). 

Self-association of SL2 

Previous works suggested that SL2 could be responsible for 
the dimerization of the viral RNA (Ivanyi-Nagy et al. 2006; 
Shetty et al. 2010). The results obtained with SL2 also 
prompted us to further investigate this stem-loop structure 
that displays a palindromic sequence. A new SL2wt was syn- 
thesized with a biotin at its 5' end. It was immobilized on a 
streptavidin-coated sensor chip. First, we checked that this 
new target could recognize a version of 5BSL3.2 shortened 
above its internal loop (Fig. IB). The result is reported in 
Figure 6. Then, SL2wt and other SL2 nonbiotinylated hair- 



pins were injected over the immobilized stem-loop SL2wt 
(Fig. 6). The results clearly show that SL2wt is able to interact 
with itself In contrast, SL2kistem) designed to generate a stable 
conformation with an apical loop exhibiting the recognition 
sequence for kissing 5BSL3.2, was not able to bind to biotiny- 
lated SL2wT- Replacement of the GU pair next to the bulge of 
conformation SL2wti (Fig- IB) by a GC pair (SL2gc) led to 
decreased binding, suggesting that increased stability of the 
3-nt stem was unfavorable for self- association. This result 
was confirmed when SL2gc was loaded on a native gel (Fig. 
3). The migration pattern of SL2gc resembles that of SL2wT) 
indicating that this stem-loop structure can still dimerize, 
but the lower bands are favored suggesting that the confor- 
mational equilibrium is shifted toward monomers. SL2gc 
behaved kineticaUy as SL2wt does for binding to 5BSL3.2: 
the sensorgrams were not monophasic (Supplemental Fig. 
S6A), which is consistent with the structural heterogeneity of 
SL2gc- Binding was abolished when, in addition, the CUA 
bulge was mutated to AAU (Fig. 6). In this case, one single 
band was observed on the native gel (Fig. 3), demonstrating 
that the bulge of SL2 was crucial for dimerization. Attempts 
to fit these data to a 1:1 model were not successful, likely due 
again to the structural heterogeneity of the SL2wt RNA im- 
mobilized on the sensor chip surface (Supplemental Fig. 
S6B-D). 

DISCUSSION 

We investigated RNA-RNA interactions at the 3' end of the 
HCV genome by SPR. This technique is a valuable tool for 
investigating molecular interactions in real time at a solid- 
liquid interface. However, it can be challenging when work- 
ing with molecules such as RNA that are prone to hydrolysis, 
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FIGURES. Kinetic analysis by the SCKODS method of SLlyst^m hairpin 
and Seq9110 binding to 5BSL3.2. The sensor chip was prepared as 
described in Figure 2 and in Materials and Methods. Seq9110 was first 
injected over the 5BSL3.2-coated surface at 1 i^M (first arrow from the 
left). During the dissociation phase of the formed Seq9110-5BSL3.2 
complex, SL2]5istem was injected sequentially in the order of increasing 
concentrations, 62.5 nM (second arrow), 250 nM (third arrow), and 
1000 nM (fourth arrow). Experiments were performed with injections 
of Seq9110 only (referred to as Seq9110 alone). The regeneration 
of the surface was achieved with a 2-min pulse of a mixture of 40% 
formamide, 3.6 M urea, and 30 mM EDTA prepared in miUi-Q water. 
The gray curves represent the experimental data (two overlaid sensor- 
grams for SL2ys,(.n, and Seq9110 alone). The black line for Seq9110 
alone represents the fit of the dissociation phase according to Equation 
(2). The black lines for SL2ystem represent the fit to a Langmuir 1:1 mod- 
el of interaction according to Equations (1) and (2) by the SCKODS 
method. 

in particular, in buffer containing magnesium ions. We have 
developed protocols that allow users to investigate complexes 
with short RNA sequences without encountering any major 
problems regarding the integrity of the RNA immobilized 
on the sensor chip surface (Di Primo et al. 2011). 

Despite previous data showing the central role of 5BSL3.2 
in the life cycle of the virus by long-range RNA-RNA inter- 
action with other genomic sequences such as SL2 (Friebe et al. 
2005; Murayama et al. 2010; Schmitt et al. 2011) and the re- 
gion centered on nucleotide 9110 (Diviney et al. 2008), sur- 
prisingly, there were no data obtained by kinetic methods 
demonstrating direct interaction between these RNA motifs. 
NMR spectroscopy has been used to characterize SL2 and 
5BSL3.2 (Friebe et al. 2005), but the kissing complex poten- 
tially formed between these two imperfect hairpins could not 
be observed. Surface plasmon resonance-based instruments 
provide real-time monitoring of virtually aU kinds of com- 
plexes since the signal results from mass changes upon bind- 
ing of an injected partner over a target-functionalized sensor 
chip surface. An interesting feature of some of these instru- 
ments resides in that the experiments can be performed at 
any temperature between 4°C and 45°C. This can be used ei- 
ther to stabilize or to destabilize complexes that could be dif- 
ficult to characterize at room temperature (23°C). During the 



preliminary experiments with SL2, we noticed that the rate of 
dissociation was very fast. We then decided to monitor the in- 
teractions at a temperature low enough (10°C) to slow down 
the dissociation phase. 

Our results clearly demonstrate that SL2 recognizes 
5BSL3.2 through kissing interactions. Interestingly, the cor- 
responding sensorgrams showed nonmonophasic shapes. 
Adding one more step to the initial Langmuir 1 : 1 model in- 
teraction would have easily fitted the data because two expo- 
nentials mathematically fit better than one, three better than 
two, etc. A typical model describing a two-step reaction is 
based on a conformational change after binding. When the 
1:1 model seems not to be appropriate, the two-step reac- 
tion model always fits the data better. The use of this model 
is subject to controversy and debate. The dogma is that SPR 
cannot detect conformational changes because the signal is 
due only to changes of mass at the solid-liquid interface. 
However, a few published works argued that only a confor- 
mational change after binding could explain unexpected sig- 
nals obtained (Winzor 2003; Wood and Lee 2005; Christopeit 
et al. 2009; Dell'Orco et al. 2010). Before we considered using 
such a model of interaction, we hypothesized that the kinetic 
characteristics of SL2wt binding to 5BSL3.2 were due to 
structural heterogeneity of the RNA sample. This was clearly 
demonstrated after loading SL2wt on a native acrylamide gel. 
Noteworthy, we could stabilize one of the predicted confor- 
mations, SL2wT2 (Fig- IB), by pairing the unpaired U with an 
A and changing the UG pair next to the apical loop to a CG 
one. The SL2wt2 conformation corresponds to that observed 
in solution by NMR spectroscopy (Friebe et al. 2005). The 
perfect hairpin that was generated from this conformation 
did follow a 1:1 model of interaction with 5BSL3.2. The exact 
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FIGURE 6. Kinetic analysis of SL2 binding to itself SL2wt was biotiny- 
lated at its 5' end and immobilized on a streptavidin-coated CMS sensor 
chip. SL2 hairpins were injected in duplicate at 500 nM as indicated by 
the arrows. The experiments were performed at 10°C and at a flow rate 
of 50 ^L/min. The truncated version of 5BSL3.2, mini5BSL3.2, was also 
injected over the SL2 functionalized surface. 
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origin of the multiphasic shape of the sensorgrams remains 
unclear (Fig. 2A). This would require further investigation 
to find which structures contribute to the signal and to 
what extent, or to analyze whether, in some cases, SPR is, in- 
deed, able to monitor conformational changes after a binding 
step. In any case, the results point out that caution should be 
taken before analyzing SPR data with kinetic models of more 
than a one-step reaction when complexity is just due to struc- 
tural heterogeneity of the samples. This is obviously of rele- 
vance when working with RNAs. 

Our data confirm previous hypotheses that suggested a 
kissing complex between 5BSL3.2 and SL2, but they do not 
agree with the results from NMR analysis (Friebe et al. 
2005) or SFiAPE mapping on the FiCV genotype lb (isolate 
Conl) (Tuplin et al. 2012). No interaction was detected by 
NMR assays with the same short RNA sequences that we 
used, and a unique secondary structure was proposed for 
SL2. We did not find a unique structure for the unmodified 
SL2. Our results demonstrate that SL2 fluctuates between at 
least two conformations: SL2wti and SL2wt2 (Fig- IB). In 
the SHAPE study on the Conl isolate, Tuplin et al. found 
that SL2 resembled the structure determined by NMR 
(SL2wT2) (Fig- IB) and that the SL2-5BSL3.2 loop-loop com- 
plex was largely absent. They concluded that the adoption of 
the known stem-loop SL2 conformation (SL2wt2) (Fig. IB) 
was incompatible with kissing- complex formation (Tuplin 
et al. 2012). This conclusion is somehow surprising. Previous 
works have shown that stable kissing complexes can only 
be formed through Watson-Crick base-pairing between the 
partially or totally complementary sequences of apical loops, 
i.e., in a stem-loop context (Westhof et al. 1988; Lee and 
Crothers 1998; Mujeeb et al. 1998; Kim and Tinoco 2000; 
Ennifar et al. 2001; Lebars et al. 2008; Van Melckebeke 
et al. 2008). The intermolecular loop-loop helix is stacked 
between the two stem helices. We observed direct interaction 
between the unmodified SL2 and 5BSL3.2, two RNA oligo- 
nucleotides similar to those used by Friebe et al. (2005). 
Furthermore, an interaction was observed that followed a 
perfect 1 : 1 model when SL2 was mutated to generate a perfect 
hairpin that displayed the proposed kissing sequence in its 
apical loop (SL2kistem) (Fig. IB). 

It is not at all clear why the NMR and SPR experiments did 
not agree even when the same RNA sequences were investi- 
gated. The NMR assays were performed at 5°C and 25°C. 
The SPR experiments were performed at 10°C. Temperature 
is unlikely to be the reason for this discrepancy. One differ- 
ence between the SPR experiments and the NMR analysis is 
that our assays were carried out in the presence of magnesium 
ions. The role of magnesium ions in the stability and struc- 
ture of RNA is well-known. The stabilizing role of magne- 
sium ions in natural or artificial kissing complexes has been 
characterized (Eguchi and Tomizawa 1991; Gregorian and 
Crothers 1995; Duconge et al. 2000). One magnesium ion 
would bind at the center of a pocket made by the two phos- 
phate clusters at the stem-loop junctions (Chang and Tinoco 



1997; Lee and Crothers 1998; Jossinet et al. 1999; Reblova 
et al. 2003). In theory, the high concentration of RNA used 
for the NMR experiments (1.2-1.4 mM) should have com- 
pensated for the absence of magnesium, since the stability 
of bimolecular nucleic acid-nucleic acid complexes depends 
on the total concentration (Puglisi and Tinoco 1989). 

The results confirm that the sequence centered on nucleo- 
tide 91 10 can recognize the internal loop of 5BSL3.2. Binding 
can also occur with an antisense sequence of the internal loop 
of 5BSL3.2, supporting previous findings that suggested that 
the sequence centered on nucleotide 9110 was likely unstruc- 
tured (Diviney et al. 2008). Together, these results support 
a view in which recognition between the internal loop of 
5BSL3.2 and the sequence centered on nucleotide 9110 
would not require structured motifs to occur, in contrast to 
what has been proposed recently (Shetty et al. 2013). The re- 
sults raise the problem of the kinetics of formation of these 
complexes. We do not understand why data obtained with 
the same immobilized target but with different partners can- 
not be fitted with a 1 : 1 model of interaction. This problem is 
further discussed in the Supplemental Material (Supplemen- 
tal Fig. S7). 

Our work demonstrates that SL2 and Seq9110 motifs do 
not compete for binding to 5BSL3.2. SL2 behaved similarly 
in the presence or in the absence of Seq9110, suggesting 
that the apical and internal loops are two structurally indepen- 
dent modules for interacting. These results, obtained with 
short RNA sequences, are in good agreement with a recent 
work with longer RNAs analyzed by SHAPE mapping that 
suggested that these interactions could occur simultaneously 
(Tuplin et al. 2012). Together, they strongly support previous 
works suggesting that 5BSL3.2 could be involved in a complex 
pseudoknot, implicating nonexclusive, long-range interac- 
tions between the internal loop of 5BSL3.2 and the sequence 
centered on nucleotide 9110 and between the apical loop of 
5BSL3.2 and SL2 (Diviney et al. 2008; Shetty et al. 2013). 

The structural heterogeneity of SL2 observed on native 
acrylamide gels and the self- association measured by SPR 
strongly suggests that SL2 can dimerize. Dimerization of viral 
genomes is normally a unique property of retroviruses in- 
volving interaction with the nucleocapsid and has been exten- 
sively documented (Paillart et al. 1996, 2004; Lu et al. 2011). 
As far as HCV is concerned, previous works have shown that 
SL2 displays a palindromic sequence, named dimer linkage 
sequence (DLS) (Cristofari et al. 2004; Ivanyi-Nagy et al. 
2006), which spans from the 5'-terminal uracile to the first 
adenine nucleotide in the apical loop of the SL2wti con- 
formation (Fig. IB). Present in the 3'X between stem-loop 
SLl and SL3, this sequence would mediate the dimerization 
of the 3' UTR in the presence of the core protein of HCV 
(Cristofari et al. 2004; Ivanyi-Nagy et al. 2006; Sharma et al. 
2012). In another work, the authors showed that the RNA 
dimerization would be mediated via a kissing complex inter- 
mediate between two structurally rearranged DLS. The core 
protein would then stabilize the resulting extended duplex 
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(Shetty et al. 2010). Our results show that the minimal struc- 
tured motif of the 3' UTR containing the dimer linkage se- 
quence SL2 is able to dimerize in vitro in the absence of 
any added protein. The results obtained with SL2]jistem and 
SL2aau bring strong experimental evidence that the upper 
bands observed in the native gel are, indeed, dimers (Fig. 
3). Upper bands are not observed if the sequence that triggers 
the dimerization is locked in a perfect stem-loop structure 
(SL2]jistem) that is too stable to allow structural rearrange- 
ments or if this sequence is mutated (SL2aau)- The biological 
role of the dimerization process is still subject to debate. It is 
commonly believed that the HCV virion contains a single 
copy of the genomic RNA, but to our knowledge, this has 
not as yet been formally established. 

We have investigated by SPR the interactions between RNA 
motifs at the 3' end of the HCV genotype lb (isolate Conl). 
We conclude that the kissing complex formed by SL2 and 
5BSL3.2, which was not observed by NMR analysis with the 
same short RNA sequences that we used or was largely absent 
when investigated by SHAPE mapping with longer RNAs, is 
compatible with the SL2wt2 conformation of SL2. This bimo- 
lecular interaction does not prevent binding of the sequence 
centered on nucleotide 9110. The two binding sequences, 
the apical and internal loops of 5BSL3.2, are structurally inde- 
pendent. The two interactions can coexist. We show that SL2 
exists in a conformational equilibrium between at least two 
conformations: one that recognizes 5BSL3.2 through loop- 
loop interactions (SL2wt2) (Fig- IB) and the other one that 
allows self-association (SL2wti) (Fig- IB). The SL2 region 
is, therefore, involved in two different interactions that are 
likely exclusive, leading consequently to different activities, 
translation, replication, or encapsidation during the viral rep- 
lication, as also suggested by a work published recently dur- 
ing the review of our manuscript (Shetty et al. 2013). 

MATERIALS AND METHODS 
Oligonucleotides 

RNA oligonucleotides, including those with a biotin tag at the 3' or 
5' ends, were synthesized on an Expedite 8908 synthesizer (Applied 
Biosystems) or were purchased from Dharmacon (Fermentas). After 
deprotection following the manufacturer's instructions, the oli- 
gonucleotides were purified by electrophoresis on 7 M urea 20% 
polyacrylamide gels, electro-eluted, and desalted on Sephadex 
G-25 (GE Healthcare Lifesciences) spin columns prepared in 1- 
mL syringes. Pure samples in miUi-Q water were assayed at 260 
nm with a Nanodrop ND-1000 spectrophotometer (Labtech) to 
determine their concentration using extinction coefficients calcu- 
lated at http://www.eurofinsdna.com/home.html. The samples 
were stored at — 20°C. Before any experiment, the RNA oligogucleo- 
tides were diluted in the appropriate buffer and heated for 1 min, 
30 sec at 90°C, then put on ice for 5 min. The predictions of the 
secondary structures were performed with the mfold web server at 
http://mfold.rna.albany.edu/?q=mfold/RNA-Folding-Form (Zuker 
2003). 



Native gels and staining of the RNA oligonucleotides 

Native gels were prepared using 15% (w/v) 75:1 acrylamide/bis(ac- 
rylamide) in 50 mM Tris-acetate buffer, pH 7.3 at 20°C, containing 
3 mM magnesium acetate (migration buffer). Two or 4 pg of RNA 
were diluted in the migration buffer, heated 1 min, 30 sec at 90°C, 
then put on ice for 5 min. The samples were loaded on the gel equil- 
ibrated at 4°C. The electrophoresis was carried out at 300 V for 5 
h. The gels were stained for 30 min in the dark using the "Stains- 
all" dye (Acros Organics), then washed with distilled water and dis- 
colored in daylight. 

Surface plasmon resonance experiments 

The SPR experiments were performed at 10°C with a Biacore T200 
apparatus (Biacore) using either SAHC200m sensor chips function- 
alized with streptavidin by the manufacturer (Xantec Bioanalytics) 
or CM5 sensor chips (Biacore). These surfaces were coated at the 
laboratory with 1000-1500 resonance units (RU) of streptavidin 
(Roche Applied Sciences) using the Biacore amine-coupling kit 
and HBS-EP running buffer (Biacore). Thirty to 70 RU of biotiny- 
lated RNA were immobilized on one flow cell of the sensor chip by 
injecting solutions prepared at 5-50 nM either in the HBS-EP buffer 
or in the running buffer used for the binding experiments (10 mM 
sodium phosphate buffer, pH 7.2 at 20°C, containing 50 mM 
sodium chloride, 3 mM magnesium chloride and 0.05% Tween- 
20). One flow cell left blank was used as a reference. All injected 
RNA samples were prepared in the running buffer. They were inject- 
ed over the functionalized surface at 50 pL/min. The binding kinet- 
ics were monitored using the single cycle kinetics (SCK) method 
(Karlsson et al. 2006), in which the samples are injected sequentially 
in the order of increasing concentrations and the regeneration step is 
performed at the end of the cycle, by the single cycle kinetics on a 
decaying surface method (Patau and Di Primo 2012), in which the 
samples are injected sequentially in the order of increasing concen- 
trations during a dissociation phase and the regeneration step is per- 
formed at the end of each cycle, or by the classical method in which 
the samples are injected one by one, with a regeneration step per- 
formed after each injection. The regeneration of the functionalized 
surface, the purpose of which is to remove any partner still interact- 
ing with the immobilized target after the dissociation phase, was 
achieved with a 2-min pulse of a mixture of 40% formamide, 3.6 
M urea, and 30 mM EDTA prepared in mUli-Q water. These condi- 
tions were mild enough so that the same functionalized sensor chip 
could be used for more than 50 cycles of injection/regeneration with- 
out major loss of the SPR signal. The sensorgrams, which correspond 
to the variation of the SPR signal expressed in resonance units (RU) 
as a function of time, were double-referenced (Myszka 1999) using 
BiaEval 4.1 (Biacore). This removes any contribution to the signal 
that may result from the RNA samples and the buffer flowing over 
the blank flow cell and the functionalized one, respectively. 

The association and dissociation rate constant and fcj, respec- 
tively, for the RNA-RNA complexes were determined by direct 
curve fitting of the sensorgrams assuming a Langmuir 1:1 model 
(see Supplemental Material for details). The dissociation equUibri- 
um constant, iCo, was calculated as fcj/fca. 

SUPPLEMENTAL MATERIAL 

Supplemental material is available for this article. 
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